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I .  Introduction 


Seals  develop  fluid-structure-interaction  forces  due  to  the 
relative  motion  of  the  rotor  and  housing  at  seal  locations/  and 
these  forces  have  been  shown  to  have  a  significant  influence  on 
the  rotordynamic  stability  and  response  characteristics  of  high 
performance  turbomachinery.  The  present  research  effort  was 
initiated  in  September  of  1982,  and  is  a  combined  computational- 
experimental  program  to  develop  experimentally-validated  compu¬ 
tational  approaches  for  predicting  seal  forces.  More  specifically 
computational  approaches  are  being  developed  to  define  the  force 
coefficients  for  the  force-motion  relationship: 


(1) 


where  (X,Y)  are  components  of  the  relative  displacement  vector 
between  the  rotor  and  the  stator  at  the  seal  location,  and 
(Fx,Fy)  are  components  of  the  reaction  force  resulting  from 
transient  pressure  and  velocity  distributions  within  the  seal. 
The  computational  approaches  which  are  being  developed  in  this 
study  will  define  the  pressure  and  velocity  distributions  within 
the  seal  which  result  from  relative  rotor-stator  motion.  The 
computational  approaches  are  to  be  validated  by  comparison  to 
experimental  measurements  of  the  following  quantities: 

(a)  transient  pressure  measurements  at  the  seal-stator 
wall  , 


(b)  transient  measurements  of  the  displacement  (X,Y)  and 


force  (Fx,Fy)  components,  which  results  in  calculated 
values  for  the  rotordynamic  (stiffness  and  damping) 
coefficients  of  Eq.  (1) ,  and 
(c)  laser  velocimeter  measurements  of  the  velocity 
distributions  within  labyrinth  seal  cavities. 

Dr.  Childs  and  Dr.  Rhode  are  responsible,  respectively,  for 
the  experimental  and  computational  developments.  In  the  following 
two  chapters,  reviews  are  provided  of  the  research  accomplishments 
with  respect  to  these  two  complementary  facets  of  the  research 
project. 

II.  Experimental  Developments 
A.  Introduction 

The  experimental  developments  of  this  program  are  jointly 
supported  by  NASA  and  AFOSR.  NASA  funds  ($223,749;  1  June  1981 
to  13  December  1983)  have  been  used  to  develop  the  test  facility 
and  test  apparatus  and  to  support  analytical-computational 
developments.  NASA  officials  are  primarily  interested  in  annular 
seals  having  constant-clearance  or  convergent-taper  geometries 
with  smooth  and  honeycomb  stator  surfaces.  Their  interest  in 
this  type  of  seal  arises  from  applications  in  the  turbine  areas 
of  the  high-pressure  turbopumps  of  the  Space  Shuttle  Main  Engine. 
This  status  report  reviews  the  experimental  progress  of  the 
complete  program  without  distinguishing  between  sources  of 
support. 


B.  Test  Apparatus 

The  test  apparatus  was  designed  and  its  manufacture  was 
supervised  by  Dr.  John  Dressman  of  the  University  of  Louisville 
who  served  as  a  subcontractor.  It  was  delivered  to  TAMU  in 
January  1983,  and  has  since  been  modified  in  the  following 
respects: 

(a)  The  housing  has  been  modified  to  accept  "swirl-rings" 
at  the  inlet.  These  swirl-rings  are  provided  to 
pre-rotate  the  entering  fluid  in  the  circumferential 
direction.  The  fluid  may  be  rotated  in  either  the 
same  direction  or  the  opposite  direction  of  rotation 
for  the  seal  rotor.  Alternatively,  straightening  vanes 
are  provided  to  minimize  fluid  prerotation.  Prior  test 
results  have  shown  that  the  inlet  tangential  velocity  is 
a  critical  factor  in  defining  the  "stiffness-cross¬ 
coupling"  coefficient  k  of  Eq.  (1) . 

(b)  The  housing  has  been  modified  to  provide  access  for 
additional  pressure  instrumentation  at  the  seal.  This 
instrumentation  includes  both  an  inlet  pitot  tube  and 
pressure  taps  along  the  seal. 

(c)  The  test  rotor  has  been  modified  to  accept  two  couplings 
which  will  allow  it  to  be  driven  by  either  an  electric 
motor  (to  5300  rpm)  or  by  an  air  turbine  (to  20,000  rpm) . 

(d)  An  additional  port  has  been  machined  in  the  housing  for 
access  by  the  LDA  system. 


(•)  Additional  holes  have  been  provided  in  the  housing  to 

provide  access  to  the  stator  for  mounting  accelerometers. 
If)  A  new  stator  has  been  manufactured  with  static  pressure 
taps  along  the  seal  and  an  inlet  pitot-tube  measurement 
for  inlet  tangential-velocity  measurement. 

(g)  A  new  high-speed  rotor  has  been  manufactured  which  can 
run  out  to  20,000  rpm. 

The  test  apparatus  utilizes  an  external  shaker  to  excite  the 
rotating  seal  within  its  stator.  The  shaker  was  purchased  using 
AFOSR  funds  ($29,000)  and  has  been  received  and  checked  out.  Some 
problems  were  encountered  at  the  outset;  however,  the  unit's  manu¬ 
facturer  (Zonic  Corporation)  has  resolved  the  problem  satisfactorily 

All  auxiliary  instrumentation,  drive-motors,  bearing-lube 
systems,  etc.,  have  been  purchased  and  checked  out. 

C.  Test  Facility 

The  test  facility  uses  an  Ingersoll-Rand  screw  compressor 
(1200  SCFM  0  100  psi)  to  supply  air  to  the  seal-test  apparatus  and 
a  Sullair  screw  compressor  (600  SCFM  @  120  psi)  to  an  air  turbine 
which  has  the  capacity  to  rotate  the  seal  up  to  20,000  rpm.  The 
compressors  were  purchased  by  TAMU  using  state  funds  and  industrial 
contributions.  Both  the  compressors  and  the  connecting  piping  to 
the  test  section  have  been  installed.  All  of  the  requisite  instru¬ 
mentations,  controls  for  the  flow  circuit,  plumbing,  electrical 
wiring,  controls,  etc.,  have  been  purchased,  installed,  and  checked 
out.  Pending  some  future  unpleasant  surprise,  which  would 
necessitate  additional  work,  the  test  facility  is  now  completed. 


D.  Graduate  Student  Development 

Mr.  Joseph  Scharrer  has  been  working  on  this  project  since 
1  January  1983.  Mr.  Scharrer  concluded  his  B.S.  degree  in 
Mechanical  Engineering  at  Northern  Arizona  University  in  May  of 
1981  (in  three  years) .  He  has  been  accepted  into  the  AFRAPT  pro¬ 
gram  by  TAMU  and  General  Electric  at  Lynn,  Massachusetts.  Mr. 
Scharrer  has  been  responsible  for  a  great  deal  of  the  experimental 
development  work  on  this  project,  and  has  recently  developed  an  im¬ 
proved  prediction  model  for  rotordynamic-coef ficients  of  a  labyrinth 
seal.  Mr.  Scharrer' s  developments  have  extended  Iwatsubo's 
analysis  [1]  to  choked  flow  conditions,  and  incorporated  a  signifi¬ 
cantly  improved  leakage  analysis  in  the  solution  procedure.  The 
analysis  developed  by  Scharrer  is  much  cruder  than  that  being 
developed  under  Dr.  Rhode's  direction.  However,  it  is  also  much 
simpler  and  provides  an  interim  prediction  tool  for  comparison  to 
experimental  test  results. 

E.  Summary 

During  the  past  year,  development  of  the  experimental  facility 
was  completed,  while  improvements  and  various  modifications  to  the 
test  apparatus  were  implemented.  An  improved  but  approximate 
analytical  solution  has  been  developed  for  the  prediction  of 
rotordynamic  coefficients  of  labyrinth  seals. 

III.  Computational  Developments 

The  serious  need  for  a  complete  recirculating  flow  model  for 
compressible  flow  labyrinth  seals  was  highlighted  at  a  recent 


workshop  [2]  held  at  TAMLJ  regarding  rotordynamic  instability 
problems.  Development  of  a  user-simplistic,  experimentally- 
validated,  finite  difference  computer  code  has  begun  for  predicting 
the  pressure  and  velocity  fields  and  subsequently  the  four 
rotordynamic  force  coefficients  of  Eq.  (1). 

The  problem  of  subsonic  air  flow  through  a  labyrinth  seal  is 
clearly  elliptic,  as  recirculation  zones  have  been  observed  in  flow 
visualization  experiments  [3].  Also,  it  is  three-dimensional 
because  of  rotor  eccentricity  effects  which  are  of  primary  concern. 
It  has  been  assumed  that  the  rotor  instability  is  manifested  as  a 
circular  precessional  orbit  of  frequency  u>  about  the  seal  center 
position.  The  desired  solution  will  be  obtained  from  the  steady- 
state  form  of  the  governing  flow  equations  because  a  steady- state 
flowfield  precesses  along  with  the  rotor  in  a  r-0  plane. 

Once  the  predicted  pressure  distribution  (for  a  given  orbital 
velocity  w)  acting  on  the  rotor  is  obtained,  it  will  be  employed 
to  calculate  the  resulting  F_  and  Fa  forces.  A  second  computation 

*r  U 

of  the  pressure  and  rotor  force  components  at  a  different  value  of 
k>  allows  both  the  stiffness  (K,k)  and  damping  (C,c)  coefficients 
to  be  calculated. 

A.  Literature  Surveys 

Substantial  progress  has  been  made  toward  the  computation  of 
rotordynamic  labyrinth  seal  coefficients.  First,  several  litera¬ 
ture  searches  were  conducted.  The  information  sought  was: 

(a)  near-wall  velocity  measurements  near  rotating  surfaces  to  allow 
derivation  of  realistic  wall  functions  needed  for  rotor  surfaces, 


(b)  recent  developments  for  reducing  false  diffusion,  and 

(c)  detailed  velocity  or  pressure  measurements  within  a  labyrinth 
seal  cavity.  The  overall  findings  of  these  searches  was  somehwat 
disappointing  in  that  the  state-of-the-art  in  the  latter  two  areas 
is  rather  primitive.  However,  useful  information  has  been 
obtained . 

B.  Computer  Plotting 

The  second  task  completed  constitutes  the  developme  .  two 
computer  software  packages  for  plotting  any  of  the  quantities  to 
be  predicted  with  the  finite  difference  code.  Both  of  these 
graphics  codes  read  labyrinth  seal  solution  values  from  computer 
disk  storage.  One  of  these  codes  plots  each  of  the  desired  radial 
profiles  at  the  proper  location  (to  scale)  within  the  corresponding 
labyrinth  seal  outline  enclosure.  The  second  graphics  package 
calculates  the  coordinates  of  contour  lines  (e.g.,  streamlines, 
isobars,  etc.)  and  plots  them  at  appropriate  locations  within  the 
seal  outline  enclosure.  At  present,  these  graphics  codes  plot 

/ 

two-dimensional  axisymmetric  (i.e.,  concentric-rotor)  predictions 
and  are  designed  for  easy  extension  to  allow  the  option  of  plot¬ 
ting  three-dimensional  (eccentric-rotor)  predictions  as  they 
become  available.  This  plotting  capability  employs  the  GCS 
graphics  package  and  is  very  helpful  in  developing  the  desired 
three-dimensional  finite  difference  code. 

C.  Computer  Code  Development 

Development  of  the  required  three-dimensional  finite  difference 
code  for  calculating  rotordynamic  force  coefficients  due  to  fluid 


forces  was  undertaken  in  a  logical,  stepwise  sequence  of  increasing 
complexity.  The  first  task  is  essentially  completed  except  for 
final  prediction  validation  tests.  This  task  consisted  of  develop¬ 
ing  a  code  for  the  two-dimensional  (axisymmetric)  compressible  flow 
in  a  single  cavity  of  a  concentric-rotor  labyrinth  seal.  This  work 
largely  consisted  of  extending  a  similar  code  previously  developed 
by  D.  L.  Rhode  for  incompressible  flows  in  labyrinth  seals,  to 
allow  the  prediction  of  corresponding  compressible  flows. 

The  next  section  discusses  the  basic  computational  methodology 
and  the  development  of  the  imcompressible  flow  code.  The  following 
two  sections  consider  the  extension  for  compressible  flow  and 
preliminary  prediction  results,  respectively. 

C.l  Incompressible  Flow  Model 

The  incompressible  flow  computer  code  for  a  labyrinth  cavity 
was  developed  from  an  improved  version  of  the  advanced  finite  dif¬ 
ference  code  recently  developed  by  Lilley  and  Rhode  [4] .  The 
computational  methodology  employed  is  advantageous  in  that  the 
task  of  extending  it  for  three-dimensional  flows  is  conceptually 
straightforward.  Numerous  modifications  to  this  code  were  required 
in  order  to  properly  incorporate  the  labyrinth  cavity  boundary 
conditions  and  flow  domain.  These  included  a  different  approach 
for  specifying  the  cavity  inflow  boundary  values  (discussed  later) 
and  the  geometric  and  boundary  condition  complexities  associated 
with  the  spinning  surface.  For  example,  it  was  decided  that  the 
code  would  have  the  capability  of  simulating  cavity  walls  of 
arbitrary  slope,  as  countless  cavity  configurations  have  been 
recently  employed  in  the  aerospace  industry. 


In  the  present  model,  all  important  physical  mechanisms  are 
simulated  by  the  complete  governing,  elliptic,  partial  differential 
equations.  These  are  the  incompressible,  axisymmetric  form  of  the 
Reynolds-averaged  equations  for  conservation  of  mass  and  momentum 
(with  x,r,6  time-mean  velocity  components  u,v,w)  as  well  as  appro¬ 
priate  turbulence  transport  equations  required  for  the  turbulence 
model.  Each  of  these  differential  equations  may  be  expressed  in 
the  general  form 

?  ♦>+  £  (Pvr  »>  -  £  <rrt  ■§£>  -  £  (rr+  ||)1  =  (2) 

where  <p  represents  any  dependent  variable  (i.e.,  u,v,w,k,  and  e), 
represents  the  appropriate  diffusion  flux  coefficient,  and  the 
source  contains  any  remaining  terms.  Symbols  k  and  e  represent 
turbulence  kinetic  energy  and  turbulence  energy  dissipation  rate. 
The  various  source  terms  for  each  equation  have  been  recently 
presented  by  Lilley  and  Rhode  and  these  are  given  in  Table  1 .  The 
standard  form  of  the  k-e  turbulence  model  as  described  by  Launder 
and  Spalding  [5]  was  employed  to  evaluate  an  assumed  isotropic 
effective  viscosity  which  is  given  by 

ueff  =  Vk2/£  +  V 

The  finite  difference  equations  are  solved  using  the  well 
known  SIMPLE  approach  [6]  on  a  complex  mesh  illustrated  in  Fig.  1. 
The  intersections,  point  P  for  example,  of  the  solid  lines  mark 
the  grid  nodes  where  all  variables  except  the  u  and  v  velocity 
components  are  stored.  The  latter  are  stored  at  points  which  are 


denoted  by  arrows  (and  labeled  w  and  s,  respectively)  located 
midway  between  the  grid  intersections,  and  the  boomerang-shaped 
envelope  encloses  a  triad  of  points  with  reference  location  P  at 
(I,J).  The  different  control  volumes  C,  U  and  V  which  are 
appropriate  for  flow  variables  stored  at  the  P,  w  and  s  locations, 
respectively,  are  given  in  Fig.  2. 

The  computational  domain  for  a  labyrinth  cavity  of  interest 
is  shown  with  a  coarse  mesh  in  Fig.  3.  There  are  several  methods 
of  numerically  simulating  the  wall  portion  which  does  not  lie 
parallel  to  a  grid  line.  In  the  interest  of  conceptual  simplicity 
the  stairstep  simulation  approach  as  previously  employed  by  Rhode, 
et  al.  [7]  was  utilized.  In  that  study  it  was  shown  that,  if  one 
is  careful,  the  error  arising  from  this  approximation  is  small  in 
comparison  to  that  of  the  turbulence  model  and  of  false  diffusion. 

The  finite  difference  equations  for  each  <j>  were  obtained  by 
integrating  Eq.  2  over  the  appropriate  control  volume  (centered 
about  the  location  of  <p)  and  expressing  the  result  in  terms  of 
neighboring  grid  point  values.  The  convection  and  diffusion  terms 
became  surface  integrals  of  the  convective  and  diffusive  fluxes 
while  the  source  term  was  linearized.  This  results  in 

'  r*  M’eAe  -  “’u*  -  r*  H]wAw 
+  [pv*  -  §*]nAn  -  tpv*  -  |4j8as  =  [sp*p  ♦  Syjvoi 

where  subscripts  n,  s,  e  and  w  refer  to  north,  south,  east  and 
west  cell  faces.  The  finite  difference  equation  for  each  <j>  cell 


became 
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Velocities  u*  and  v*  and  pressure  p*  are  only  estimates  prior 
to  convergence  and,  for  example,  the  axial  momentum  equation  which 
is  solved  for  u*  is 


u  *  u  *  .  * 
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The  resulting  local  value  of  u*  must  be  corrected  in  order  to  also 
be  an  approximate  local  solution  of  the  mass  conservation  equation. 
This  correction  occurs  via  the  pressure  correction  P '  (calculated 
from  mass  and  momentum  conservation  considerations)  according  to 


u 


=  u 


dU<pw 


"  PP> 


where  Du  =  i^/ap.  The  radial  velocity  v  is  similarly  corrected, 
whereas  pressure  is  corrected  by 


The  P'  equation  is  a  type  of  Poisson  pressure  equation  given  by 


appp  =  l  ajpj  +  sl 


Numerical  solution  proceeds  by  iterating  the  following  steps: 

(1)  guess  values  of  all  variables;  (2)  calculate  u*  and  v*  from  the 
momentum  difference  equations;  (3)  solve  the  pressure  correction 
equation  for  P";  (4)  obtain  pressure  P  and  corrected  velocities  u 
and  v;  (5)  solve  other  difference  equations  for  other  <j>  variables 
successively;  and  (6)  return  to  (2)  using  the  new  values  as 
improved  guesses. 

The  numerical  scheme  considers  variables  at  points  along  a 
vertical  grid  line  as  unknowns  (e.g.,  N,  P,  and  S  for  each  point  P) 
whereas  the  latest  values  of  each  E  and  W  neighbor  are  treated  as 
known.  The  tridiagonal  matrix  algorithm  (TDMA)  is  employed  to 
sequentially  solve  the  system  of  algebraic  equations  resulting  for 
each  vertical  line,  sweeping  from  left  to  right. 

The  boundary  conditions  were  numerically  incorporated  via 
appropriate  modification  of  the  finite  difference  equation  for 
computational  cells  adjacent  to  a  boundary.  Boundary  values  at 
the  seal  cavity  flow  inlet  are  naturally  very  important,  but 
unfortunately  remain  unknown,  as  necessary  quantitative  flow 
measurements  do  not  exist  even  at  a  single  operating  condition. 
Current  plans  here  at  Texas  A&M  University  include  extensive 
measurement  of  the  necessary  values  using  an  advanced  laser  dop- 
pler  anemometer.  Due  to  this  current  lack  of  inlet  boundary 
data,  current  practice  is  that  for  each  iteration  the  inlet 


values  of  each  variable  (except  pressure)  are  assumed  equal  to 
the  latest  corresponding  outlet  values.  This  practice  implies 
a  series  of  geometrically  identical  cavities  (with  a  non-stepped 
stator  housing)  and  a  streamwise  periodic  flow  which  occurs 
downstream  of  the  first  two  or  three  cavities. 

No  boundary  conditions  are  needed  at  the  cavity  outlet  be¬ 
cause  of  the  upwind  differencing  formulation  which  occurs  there 
for  all  variables.  The  x-  and  r-momentum  equations  require  the 
zero  convective  and  diffusive  flux  specification  along  walls  which 
are  normal  to  the  respective  velocity  component  directions. 

Boundary  values  for  these  two  equations  along  walls  which  are 
parallel  to  respective  velocity  component  directions  are  implied 
through  the  use  of  experimentally-based  wall  functions.  Appropriate 
wall  functions  are  similarly  employed  for  the  0-momentum  equation 
along  all  impervious  walls.  Wall  functions  are  also  used  along 
wall  boundaries  for  the  k-  and  e-equations  as  derived  for  swirling 
flows  by  Lilley  and  Rhode  [4], 

C.2  Compressible  Flow  Model 

The  incompressible  flow  code  described  above  was  extended 
to  allow  subsonic  compressible  flow  computations  in  a  labyrinth 
cavity.  Four  primary  modifications  were  implemented  to  accomplish 
this.  First,  the  fluid  was  assumed  to  be  an  ideal  gas  (air)  and 
an  expression  for  calculating  density  (the  ideal  gas  law)  was 
incorporated. 

The  second  modification  consisted  of  including  another  simul¬ 
taneous  partial  differential  equation  (conservation  of  energy) 


from  which  the  temperature  T  would  be  determined.  The  energy 
equation  in  terms  of  stagnation  enthalpy  H  is 


757  =  dr  ‘v-j  +  K  dr» 


where 


H  =  Cp  T  +  \  V 


and  the  expression  for  is  given  subsequently.  This  equation 
was  expressed  as 
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wherein  Cp  was  assumed  constant. 

The  general  form  of  this  equation  which  was  incorporated  is 
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The  enthalpy  boundary  conditions  included  the  specification 
of  inlet  temperature  values  which  must  be  known  or  estimated,  and 
all  walls  were  treated  as  adiabatic.  Boundary  conditions  for  the 
kinetic  energy  portion  of  H  were  specified  through  wall  functions 
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as  for  the  momentum  equations. 

The  stres  tensor  was  augmented  to  account  for  compres¬ 
sibility  effects  in  the  usual  way  as  the  third  modification.  The 
resulting  stress  relation  is  given  by 


Tij 


^eff  (9vi/3xj  +  3V3V  "  6ij  3(pk  +  pef f  9vm/axm) 


where  the  underlined  quantity  is  the  new  contribution,  as  it  is 
not  generally  negligible  for  compressible  flows.  The  alteration 
was  effected  where  necessary,  i.e.,  in  the  finite  difference 
equations  for  u,  v,  k,  e  and  H. 

The  fourth  modification  involved  the  derivation  and  imple¬ 
mentation  of  the  compressible  flow  form  of  the  equation  for 
pressure  correction  P';  a  literature  search  had  revealed  that 
this  equation  is  not  available.  Some  highlights  of  the  derivation 
are  briefly  considered  here  for  the  simple  case  of  one-dimensional 
social  flow.  The  derivation  began  with  the  discretized  form  of  the 
steady  flow  mass  conservation  equation  which  is 


(pu)e  -  (pu)w  =  0 


(3) 


where  the  lower  case  subscripts  refer  to  values  at  control  volume 
faces. 

A  typical  control  volume  for  P'  is  centered  at  a  typical  P' 
storage  location,  i.e.,  a  point  labeled  P  in  Figs,  land  2.  It 
is  the  C  control  volume  which  is  utilized  for  the  P'  equation. 

The  east  face  mass  flux,  for  example,  was  expressed  as 


(pu)e  =  (p*+p')e  (U*+U')e  ^ 


'V  p*u*  +  p*u' 


e  e 


+  p'u* 
e  e  e  e 


where  the  term  p^u^  was  neglected.  It  was  assumed  that 

pe  =  J5<pe+pP>  =  ^(pe/rtE  +  pp/RTp) 

using  the  ideal  gas  law  (where  capitalized  subscripts  refer  to 
grid  point  values).  From  axial  momentum  considerations,  one  finds 
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as  given  earlier. 

Upon  substituting  the  appropriate  east  and  west  face  mass 
flux  expressions  into  Eq.  (3),  the  following  was  obtained  for 
the  P'  difference  equation 

{AE  +  AW  +  F]  Pp  =  AE  Pjj,  +  AW  PT^  +  Am* 


where 
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AE  =  (p*  -  u|/2RTe) 

AW  =  (p*  d"  +  u*/2RTw) 

F  =  (u*/2RT£  +  U*/2RTp)  -  (u*/2RTw  +  u*/2RTp) 

Am*  =  p*u*  -  p*u* 

The  P'  equation  actually  implemented  was  derived  from  the  two- 
dimensional  flow  mass  conservation  equation,  although  the  procedure 
exactly  parallels  that  given  here.  The  appropriate  boundary  con¬ 
dition  which  was  implemented  for  this  equation  is  P'=0  along  all 
boundaries. 

The  predictive  capability  of  the  compressible  flow  code  remains 
to  be  evaluated  for  a  compressible  labyrinth  seal  flow.  This 
results  from  the  lack  of  appropriate  measurements  for  comparison. 
Measurements  of  well-documented  similar  flows  are  currently  being 
sought  for  this  purpose.  However,  a  small  sample  of  incompressible' 
flow  labyrinth  seal  measurements  were  found;  these  are  employed 
in  the  next  section  in  a  preliminary  prediction  test  case.  The 
next  section  also  includes  detailed  predicted  results  of  another 
incompressible  flow  in  a  labyrinth  seal,  as  these  are  of  considerable 
interest. 


N 


D.  Preliminary  Predictions  of  the  Concentric-Rotor  Seal 
D.l  Incompressible  Flow  Test  Case 
A  computerized  literature  search  indicated  that  only  the 
paper  by  Stoff  [8]  gives  detailed  quantitative  velocity  measure¬ 
ments  within  a  labyrinth  seal  cavity.  Only  the  circumferential 


1 


mean  velocity  profile  is  presented;  it  was  measured  with  a  laser 
doppler  anemometer  at  a  single  axial  station  located  midway 
between  adjacent  teeth.  These  measurements  were  employed  for 
preliminary  verification  testing  of  the  current  computer  code, 
although  much  more  data  is  needed  for  final  testing.  As  previously 
mentioned,  current  plans  are  that  G.  L.  Morrison  will  obtain  the 
required  measurements  here  using  a  laser  doppler  anemometer 
purchased  with  Defense  Department  funds. 

Stoff  utilized  a  large-scale  model  test  section  with  identical 
(non-stepped)  rectangular  cavities  through  which  water  flowed. 

It  is  assumed  that  the  flow  conditions  are  characterized  primarily 
by  the  axial  and  circumferential  Reynolds  numbers  as  well  as  the 
Taylor  number.  These  are  based  on  spatially  averaged  swirl  velo¬ 
city  W,  tooth-clearance  average  axial  velocity  0,  kinematic 
viscosity  v,  shaft  radius  rgh,  tooth  clearance  c,  and  radial 
distance  from  the  shaft  to  the  stator  wall  d; 


Rex  -  T-  %  1.7  x  102 


Re0  =  58-  %  3.0  x  104 


Ta  -  T  S  ia  * io4 


Figure  4  shows  that  predictions  using  the  incompressible  flow 
computer  code  are  in  good  agreement  (6%  discrepancy)  with  Stoff's 
corresponding  measurements.  Further,  Stoff's  computation  employed 
the  same  numerical  method  and  hybrid  differencing  scheme  discussed 
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earlier  and  his  predictions  are  essentially  identical  to  those 
shown  here. 

D.2  Preliminary  Predictions 

The  predicted  results  for  the  concentric-rotor  case  of  an 
incompressible  labyrinth  seal  flow  field  are  included  here.  This 
serves  to  enhance  one's  understanding  of  the  eccentric-rotor,  com¬ 
pressible  flow  case  which  produces  rotor  instability  forces.  It 
is  the  "interstage"  seal  of  the  liquid  hydrogen  turbopump  for  the 
Space  Shuttle  main  engine  which  is  considered.  Predictions  were 
obtained  for  both  a  generic  as  well  as  a  current  seal  design.  Both 
of  these  designs  exhibit  a  non-stepped  stator  wall.  As  noted 
earlier,  this  design  feature  alleviates  the  difficulty  that  the 
inlet  flow  boundary  conditions  are  not  known;  i.e.,  it  allows  one 
to  specify  inlet  values  of  each  variable  (except  pressure)  equal 
to  corresponding  outlet  values  for  each  numerical  iteration.  This 
results  from  the  fact  that  the  cavity  inlet  and  outlet  tooth 
clearance  areas  are  the  same.  This  practice  assumes  that  the 
flow  through  a  series  of  cavities  is  streamwise  periodic. 

All  operating  and  geometrical  conditions  are  typical  of  actual 
flight  conditions.  Liquid  hydrogen  at  42 °R  enters  the  cavity  with 
tooth-clearance  spatially  averaged  mean  velocity  U  =  338  m/s.  The 
primary  dimensionless  flow  parameters  (defined  earlier)  are: 

(a)  Re  =  5.0  x  105,  (b)  Ren  =  7.6  x  105,  and  (c)  Ta  =  1.2  x  105. 

X  U 

-4 

The  shaft  speed  ft  =  35,410  rpm,  cavity  axial  length  L  =  8.585  x  10  m, 

_2 

and  stator  wall  radius  R  =  4.29  x  10  m. 

A  grid  dependence  study  was  conducted  which  revealed  that  a 
uniform  53  x  53  grid  is  satisfactory,  exhibiting  negligible  grid 
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size  effects.  Convergence  tests  using  convergence  criteria  as 
low  as  0.05%  for  the  normalized  sum  of  residual  mass  source 
magnitudes ,  for  example,  were  analyzed. 

Figure  5  shows  predicted  dimensionless  streamline  patterns 
calculated  and  plotted  by  computer.  The  recirculation  cell  given 
in  Fig.  5(a)  is  the  expected  result  of  a  free  shear  layer  flowing 
over  a  cavity.  The  predictions  give  no  indication  that  an  addi¬ 
tional  small  recirculation  cell  exists  in  the  corner  at  the  base 
of  each  tooth.  Note  that  for  the  case  of  negligible  inlet  (i.e., 
leakage)  flow  rate,  a  pair  of  counter-rotating  Taylor-like  cells 
are  expected  in  the  cavity.  Stoff  presents  a  predicted  velocity 
vector  plot  of  this  case.  Also,  he  shows  a  similar  plot  for  the 
intermediate  case  of  small  inlet  flow  rate  U/W  =  0.06.  This  flow 
exhibits  both  recirculation  cells,  the  upstream  one  being  extended 
downstream  near  the  free  shear  layer,  and  the  other  occupying  only 
the  lower  right-hand  corner  of  the  rectangular  cavity. 

Also  observe  in  Fig.  5(a)  that  the  dividing  streamline  has  a 
stagnation  point  slightly  below  the  corner  of  the  downstream  tooth. 
This  is  also  true  of  the  typical  cavity  design  as  seen  in  Fig.  5(b). 
In  the  latter  case  the  recirculation  zone  is  centered  slightly 
downstream  and  above  that  for  the  generic  cavity.  Also,  the  pre¬ 
dicted  recirculation  strength  is  reduced  to  approximately  29%  of 
the  leakage  flow  rate,  whereas  it  is  39%  for  the  generic  design. 

As  previously  mentioned,  many  quantitative  details  for  the 
flow  field  in  a  labyrinth  cavity  are  not  found  elsewhere,  particu¬ 
larly  concerning  the  inlet  boundary  values.  Therefore  radia 
profiles  of  numerous  flow  variables  are  included  here.  Figui  6(a) 


shows  the  dimensionless  mean  axial  velocity  distribution  for  the 
generic  cavity.  A  very  large  value  of  3u/3r  occurs  near  the  inlet 
separation  corner;  this  is  indicat-./'i'e  of  high  turbulence  energy 
generation  in  the  free  shear  layer.  The  maximum  velocity  at  a 
given  axial  station  occurs  very  near  the  stator  wall  due  to  the 
large  centrifugal  force.  The  axial  velocity  distribution  for  the 
typical  cavity  in  Fig.  6(b)  is  very  similar.  The  most  significant 
differences  appear  at  x/L  =  0.35  near  the  free  shear  layer  and  in 
the  base  of  the  cavity  where  the  flow  along  the  shaft  is  reduced. 

As  expected,  there  is  negligible  mean  radial  velocity  in  the 
leakage  flow  region  as  seen  in  Fig.  7(a),  except  near  the  stagna¬ 
tion  point  on  the  downstream  tooth.  This  portion  of  the  leakage 
flow  is  slightly  increased  for  the  typical  cavity  of  Fig.  7(b). 

Seal  designers  can  probably  reduce  this  leakage  contribution  by 
altering  the  tooth  shape  at  this  location.  The  swirl  velocity 
in  Figs.  8(a)  and  (b)  is  normalized  by  the  shaft  circumferential 
velocity.  The  values  for  both  seal  designs  are  almost  uniform 
throughout  the  entire  domain,  with  a  sharp  increase  very  near  a 
rotating  surface. 

Figures  9(a)  and  (b)  show  the  distribution  of  pressure  relative 
to  the  cavity  inlet  stator  wall  pressure  PQW-  In  the  leakage  flow 
region  of  the  generic  case,  pressure  decreases  in  the  streamwise 
direction  until  about  x/L  =  0.5.  Then  it  increases  slightly  until 
approximately  x/L  =  0.78,  whereupon  it  suddenly  decreases  as  the 
flow  accelerates  over  the  downstream  tooth.  This  pressure  varia¬ 
tion  is  a  reflection  of  the  Bernoulli  effect  of  pressure  decrease 
at  constrictions  in  flow  area.  These  variations  of  flow  area  are 


indicated  by  the  dividing  streamline  in  Fig.  5(a).  Within  the 
recirculation  region  also,  pressure  exhibits  the  expected  behavior 
with  the  minimum  value  at  the  center  of  the  vortex.  A  sharp 
positive  pressure  peak,  which  occurs  in  an  annular  fashion  at 
x/L  =  0.85,  results  from  the  flow  stagnation  on  the  downstream 
tooth.  The  resulting  overall  pressure  drop  from  inlet  to  outlet 
of  this  cavity  at  the  stated  conditions  is  predicted  as  600  kPa. 

The  pressure  distribution  within  the  typical  cavity  is  seen 
if  Fig.  9(b)  to  exhibit  similar  trends.  However,  leakage  flow 
pressure  varies  in  the  streamwise  direction  much  more  gradually 
for  this  geometry,  until  about  x/L  =  0.65.  Proceeding  downstream, 
pressure  decreases  more  rapidly  than  for  the  generic  case  as  the 
flow  area  constricts  in  passing  over  the  downstream  tooth.  Another 
effect  of  cavity  design  on  pressure  is  the  reduced  level  of  3P/3r 
within  the  typical  cavity,  a  consequence  of  the  reduced  recircula¬ 
tion  strength  mentioned  earlier.  The  predicted  bulk  pressure  drop 
for  this  cavity  is  approximately  an  increase  of  30%  from  the 
generic  geometry.  Thus,  these  pressure  predictions  reveal  that 
the  typical  cavity  is  indeed  preferred  over  the  generic  one. 

As  mentioned  earlier,  the  large  value  of  3u/3r  in  the  free 
shear  layer  yields  a  high  generation  rate  of  turbulence  kinetic 
energy,  which  of  course  promotes  the  large  overall  pressure  drop. 
Hence,  the  station  maximum  of  turbulence  energy  occurs  near  the 
free  shear  region  as  seen  in  Figs.  10(a)  and  (b) .  Note  that  the 
typical  cavity  exhibits  a  significantly  higher  level  of  turbulence 
energy  than  for  the  generic  case.  The  greatest  value  occurs  in 
both  cases  at  x/L  =  0.85  near  the  stagnation  point.  A  significant 


portion  of  this  turbulence  energy  is  convected  up  and  over  the 
downstream  tooth  as  indicated  by  the  distribution  at  x/L  =  1.0. 

The  remainder  is  largely  dissipated  by  viscous  action  in  the 
recirculation  zone. 

The  accuracy  of  the  present  predictions  is  difficult  to 
assess  due  to  the  severe  lack  of  suitable  measurements  for  bench¬ 
mark  prediction  tests  of  seal  flows.  The  highly  successful 
prediction  of  Stoff's  swirl  velocity  measurements  does,  however, 
provide  an  indication  of  considerable  realism.  The  present  pre¬ 
dictions  do  exhibit  some  degree  of  false  diffusion  (numerical) 
error,  but  this  is  largely  limited  only  to  the  tip  and  base  region 
of  each  tooth.  An  additional  source  of  error  for  most  swirling 
flow  is  inherent  in  the  k-e  turbulence  model,  and  thus,  further 
work  is  needed  in’  developing  a  more  accurate  approach  to  turbulence 
closure. 

E .  Summary 

Progress  in  developing  the  necessary  three-dimensional  flow 
computer  code  for  predicting  the  rotordynamic  force  coefficients 
included:  (a)  several  computerized  literature  searches, 

(b)  development  of  two  graphics  packages  for  computer  plotting 
of  predictions,  (c)  development  of  the  task-one  compressible  seal 
flow  prediction  code,  (d)  preliminary  prediction  validation  test¬ 
ing  using  Stoff’s  limited  measurements,  and  (e)  brief  exploratory 
predictions  to  obtain  additional  insight  regarding  unknown  seal 
flowfield  boundary  conditions,  etc.  The  task-one  code  allows  the 
prediction  of  the  two-dimensional,  axisymmetric  air  flow  through 
a  cavity  of  a  concentric-rotor  labyrinth  seal. 


Generation  of  this  code  involved  numerous  developments 
applied  to  an  existing  two-dimensional  flow  code.  These  included: 
(a)  implementation  of  many  computational  domain  boundary  and 
boundary  condition  complexities,  (b)  incorporation  of  the  conser¬ 
vation  of  energy  differential  equation,  (c)  addition  of  many  stress 
terms  resulting  from  compressibility  effects,  ?nd  (d)  derivation 
and  proper  implementation  of  a  new  pressure  correction  equation. 

F.  Graduate  Student  Development 

Mr.  Steve  Sobolik  has  been  employed  on  this  project  since 
1  September  1982.  He  received  a  B.S.  degree  (Cum  Laude)  in 
Mechanical  Engineering  from  Texas  A&M  University  in  May  1982.  He 
designed  and  developed  software  for  Mitre  Corp.  (a  NASA  contractor) 
during  three  summers.  Mr.  Sobolik  has  been  responsible  for  most 
of  the  computational  development  work  on  this  project.  He  is  an 
innovative  individual  and  produces  quality  work. 

G.  Publications 

A  final  manuscript  concerning  this  investigation  is  currently 
in  preparation  as  a  technical  journal  article.  Current  plans  are 
that  it  will  be  submitted  to  the  ASME  Journal  of  Fluids  Engineering. 
The  tentative  title  is:  "On  the  Prediction  of  Incompressible  Flow 
in  Labyrinth  Seals”.  Tentative  authors  include  D.  L.  Rhode,  S.  R. 
Sobolik,  and  J.  A.  Demko. 

H.  Conference  Presentations 

The  above  paper  has  been  accepted  for  presentation  at  the 
ASME  Fluids  Engineering  Conference  to  be  held  in  New  Orleans, 
February  11  through  17,  1984.  Also,  an  abstract  entitled 


"Prediction  of  Compressible  Flow  in  a  Labyrinth  Seal"  has  been 

submitted  to  the  Fluid  Dynamics  Division  of  the  American  Physical 

Society.  The  meeting  will  be  held  in  Houston,  November  20  through 

22,  1983.  Tentative  authors  are  D.  L.  Rhode  and  S.  R.  Sobolik. 
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THREE  GRIDS  :  FOR  p,w  ETC.  -  AT  POSITION  MARKED  (O) 

FOR  u  VELOCITY  -  AT  POSITION  MARKED  (—) 
FOR  v  VELOCITY-  AT  POSITION  MARKED  (  ♦  ' 


Fig.  1.  Staggered  grid  and  notation  for  the  rectangular  computational 


(b)  Typical  cavity 


Predicted  dimensionless  mean  radial  velocity  for: 
(a)  a  generic  and  (b)  a  typical  labyrinth  seal  cavi 


